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Abstract. We examine the non-ergodic properties of scaled Brownian motion, a non¬ 
stationary stochastic process with a time dependent diffusivity of the form D{t) ~ 

We compute the ergodicity breaking parameter EB in the entire range of scaling 
exponents a, both analytically and via extensive computer simulations of the stochastic 
Langevin equation. We demonstrate that in the limit of long trajectory lengths T and 
short lag times A the EB parameter as function of the scaling exponent a has no 
divergence dX a = ll2 and present the asymptotes for EB in different limits. We 
generalise the analytical and simulations results for the time averaged and ergodic 
properties of scaled Brownian motion in the presence of ageing, that is, when the 
observation of the system starts only a finite time span after its initiation. The 
approach developed here for the calculation of the higher time averaged moments 
of the particle displacement can be applied to derive the ergodic properties of other 
stochastic processes such as fractional Brownian motion. 


PAGS numbers: 05.40.a,02.50.-r,87.10.Mn 


1. Introduction 


The non-Brownian scaling of the mean squared displacement (MSD) of a diffusing 
particle of the power-law form PEIEHU 

{x\t)) = 2K^r (1) 


is a hallmark of a wide range of anomalous diffusion processes Ell]. Equation (fT|) 
features the anomalous diffusion coefficient Ka of physical dimension cm^/sec" and 
the anomalous diffusion exponent a. Depending on its magnitude we distinguish 
subdiffusion (0 < a < 1) and superdiffusion {a > 1). Interest in anomalous 
diffusion processes was rekindled with the advance of modern spectroscopic methods, 
in particular, advanced single particle tracking methods [5]. Thus, subdiffusion was 


Quantifying the non-ergodicity of scaled Brownian motion 


2 


observed for the motion of biopolymers and submicron tracer particles in living biological 
cells [6], in complex fluids [7], as well as in extensive computer simulations of membranes 
[8] or structured systems [9], among others [31 lU [ID]. Superdiffusion of tracer particles 
was observed in living cells due to active motion m- 

Anomalous diffusion processes characterised by the MSD ([T]) may originate from 
a variety of distinct physical mechanisms HI m n uni ini US]. These include a power- 
law statistic of trapping times in the continuous time random walks (CTRWs) as well 
as related random energy models n [ini [El [isi n 15] and CTRW variants with 
correlated jumps [16] or superimposed environmental noise HZI. Other models include 
random processes driven by Gaussian yet power-law correlated noise such as fractional 
Brownian motion (FBM) [18] or the fractional Langevin equation [19]. Closely related 
to these models is the subdiffusive motion on fractals such as critical percolation clusters 
ra- Finally, among the popular anomalous diffusion models we mention heterogeneous 
diffusion processes with given space dependencies of the diffusion coefficient [21] as 
well as processes with explicitly time dependence diffusion coefficients, in particular, 
the scaled Brownian motion (SBM) with power-law form D{t) ~ analysed in more 
detail herein [22] [23] [2ll [25] . Also combinations of space and time dependent diffusivities 
were investigated [23] [26]. Space and/or time dependent diffusivities were used to 
model experimental results for smaller tracer proteins in living cells [27] and anomalous 
diffusion in biological tissues [28] including brain matter [29] [30]. In particular, SBM 
was used to describe fluorescence recovery after photobleaching in various settings [31] 
as well as anomalous diffusion in various biophysical contexts [32] . In other branches of 
physics SBM was used to model turbulent flows observed by Richardson [33] as early 
as 1952 by Batchelor [34]. Moreover, the diffusion of particles in granular gases with 
relative speed dependent restitution coefficients follow SBM [35]. We note that in the 
limiting case D{t) ~ 1/t the resulting process is ultraslow with a logarithmic growth 
of the MSD [36] known from processes such as Sinai diffusion [3H], single hie motion in 
ageing environments [39], or granular gas diffusion with constant restitution coefficient 

ra. 

In the following we study the ergodic properties of SBM in the Boltzmann-Khinchin 
sense [37] , hnding that even long time averages of physical observables such as the MSD 
do not converge to the corresponding ensemble average [1] [E] [E] HD]- In particular 
we compute the ergodicity breaking parameter EB—characterising the trajectory-to- 
trajectory huctuations of the time averaged MSD—in the entire range of the scaling 
exponents a, both analytically and from extensive computer simulations. We generalise 
the results for the ergodic properties of SBM in the presence of ageing, when we start to 
evaluate the time average the MSD a hnite time span after the initiation of the system. 

The paper is organised as follows. In section [2] we summarise the observables 
computed and provide a brief overview of the basic properties of SBM. In section [3] 
we describe the theoretical concepts and numerical scheme employed in the paper. We 
present the main results for the EB parameter of non-ageing and ageing SBM in detail 
in sections [3] and [4] In section [5] we summarise our Endings and discuss their possible 
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2. Observables and fundamental properties of scaled Brownian motion 


We define SBM in terms of the stochastic process [H [22l [2ll [26l 03] 

^ = Vmi) X C(i), (2) 

where ({t) is white Gaussian noise with zero mean and unit amplitude (C(G)C(G)) = 
6{ti — ^ 2 )- The time dependent diffusion coefhcient is taken as 

D{t) = aKj^-\ (3) 


where we require the positivity of the scaling exponent, a > 0. SBM is inherently out of 
thermal equilibrium in conhning external potentials 12 a. Let us briefly outline the basic 
properties of the SBM process. The ensemble averaged MSD of SBM scales anomalously 
with time in the form of equation (JI]). 

Here and below we use the standard dehnition of the time averaged MSD mm 


T-A 


(52(A) = 


T-A 


x{t + A) — x{t) 


1 2 


dt, 


(4) 


where A is the lag time, or the width of the window slid along the time series in taking 
the time average (jl]). Moreover, T is the total length of the time series. We denote 
ensemble averages by the angular brackets while time averages are indicated by the 
overline. Often, an additional average of the form 

, N 




(5) 


2=1 


is performed over N realisations of the process, to obtain smoother curves. From a 
mathematical point of view, this trajectory average allows the calculation of the time 
averaged MSD for processes, which are not self-averaging [H Both quantities (0]) 
and (EJ) are important in the analysis of single particle trajectories measured in advanced 
tracking experiments [12]. For SBM the mean time averaged MSD ([5]) grows as 

^0+1 _ ^«+i _ _ A)"'*'^] 




(a + l)(r-A) 


( 6 ) 


In the limit A/T <C 1, the time averaged MSD scales linearly with the lag time, 

~ 2AA (7) 

SBM is thus a weakly non-ergodic process in Bouchaud’s sense HD: the ensemble 
and time averaged MSDs are disparate even in the limit of long observation times T, 
limT^oo <5^(A) ^ (a;2(t)) and thus violate the Boltzmann-Khinchin ergodic hypothesis, 
while the entire phase space is accessible to any single particle. Moreover, the magnitude 


That is, a sufficiently long time average is sufficient to represent the whole ensemble. 
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of the time averaged MSD becomes a function of the trace length T. Analogous 
asymptotic forms for the mean time averaged MSD ([5]) are found in subdiffusive CTRW 
processes [iQllil] and heterogeneous diffusion processes IZB, see also the extensive recent 
review [1]. Note that also much weaker forms of non-ergodic behaviour exist for Levy 
processes m 

Another distinct feature of weakly non-ergodic processes of the subdiffusive CTRW 
[40] and heterogeneous diffusion type [21] is the fact that time averaged observables 
remain random quantities even in the long time limit and thus exhibit a distinct scatter 
of amplitudes between individual realisations for a given lag time. This irreproducibility 
due to the scatter of individual traces <5^ (A) around their mean is described by the 
ergodicity breaking parameter jH HD] SS] 146] 


EB(A) 




V(A) 

D(A) 


(«"(A)> - 1, 


( 8 ) 


where ^(A) = ^^(A) j ^(52(A)y. Moreover, we introduced the abbreviations A/'(A) and 
T>{A) for the nominator and denominator of EB, respectively. This notation will be 
used below. For Brownian motion in the limit A/T —)■ 0 the EB parameter vanishes 
linearly with A/T in the form jH 145] 


EBbm(A) - —. 


(9) 


In contrast to subdiffusive CTRW and heterogeneous diffusion processes, the EB 
parameter of SBM vanishes in the limit A/T ^ 0 and in this sense the time averaged 
observable becomes reproducible [24] |25l |43]. We demonstrate the small amplitude 
scatter of SBM in hgure [U for a detailed discussion see below. We note that the scatter 
of the time averaged MSD of SBM around the ergodic value = 1 becomes progressively 
asymmetric for smaller a values and in later parts of the time averaged trajectories, see 
Fig. 6 of reference [4]. In the following we derive the exact analytical results for the 
EB parameter of SBM and support these results with extensive computer simulations. 
Moreover we extend the analytical and computational analysis of the EB parameter to 
the case of the ageing SBM process when we start evaluating the time series x(t) at the 
time ta > 0 after the original initiation of the system at t = 0 [43] . 

The time averaged MSD of an ageing stochastic process is dehned as [T5] 


rta+T-A 


5^(A) = 


'ta 


x{t -I- A) — x{t) 


dt 


( 10 ) 


T-A 

and thus again involves the observation time T. The properties ageing SBM were 
considered recently [43]. The mean time averaged MSD becomes 


<5^(A) = 


2JC 


(T + 4)“+1-(4 + A) 


(« + l)(T-A)L 

-(T + 4-Ar+i 


a+l 


+ t' 


a+l 


(11) 
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Figure 1. Time averaged MSD of SBM as function of the lag time A for several values 
of the scaling exponents a and ageing times ta- The asymptotic behaviour of equation 
(fTT|l is shown by the black solid lines. Parameters: T = 10^, ta = 0, 10^, 10^, and 
N = 100 traces are shown. 


The ratio of the aged versus the non-ageing time averaged MSD in the limit A <^ta,T 
has the asymptotic form [13] 

(^) 

Aa(ta/T) = ) ( ~ (1 + ta/rr " (ta/T)“ (12) 

(<52(A)) 

This functional form is identical to that obtained for subdiffusive CTRWs [15] and 
heterogeneous diffusion processes [17|. The factor Aa{z) quantifies the respective 
depression and enhancement of the time averaged MSD for the cases of ageing sub- 
and superdiffusive SBM. 

Figured] shows the time averaged MSD (A) of individual SBM traces for the case 
of weak, intermediate, and strong ageing for different values of a. We observe that the 
spread of individual 5^(A) changes only marginally with progressive ageing times ta- 
Also the changes with the scaling exponent a are modest, compare hgureO Also note 
that the magnitude of the time averaged MSD decreases with ta for ultraslow SBM at 
a = 0, stays independent on ta for Brownian motion at a = 1, and increases with the 
ageing time for superdiffusive processes at a > 1. These trends are in agreement with 
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Figure 2. Distribution of the relative amplitude S‘^{A) j of the time 

averaged MSD traces for SBM processes with different scaling exponents a as indicated 
in the panels. As expected, the spread grows and the distribution becomes more 
leptokurtic at longer lag times A. For progressively larger values of the scaling 
exponent a the spread of time averaged MSD decreases but stays asymmetric with 
a longer tail at larger 6^ values. In particular, for a = 1 and 2 the shape is almost 
indistinguishable at A = 10, see the bottom right panel. The trace length is T = 10“^ 
and the number of traces used for averaging is 10^. 


the theoretical predictions of equation flTT]) shown as the solid lines in hgure [T] 

3. Ergodicity breaking of non-ageing scaled Brownian motion 
3.1. General expression for the ergodicity breaking parameter 

Analytically, the derivation of the EB parameter for SBM involves the evaluation of the 
fourth order moment of the time averaged MSD, 

T-A T-A 

^((52(A)) ^ j dti j dt 2 (^(x\ti + A) - x{G)f 

0 0 

x(x^((2 + A) - ^(( 2 )) 



(13) 
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We use the fundamental property of SBM that 

{x{ti)x{t2)) = (a:^(min{ti,f2})) , (14) 

and the Wick-Isserlis theorem for the fourth order correlators [H]. We then obtain the 
nominator J\f of the EB parameter of equation (0) 

Af(A)= 

rT-A 


-T-A 


(T - A)^ 


dti j dt 2 (^{x{ti +/^) - x{ti)) 
X {x{t 2 + A) - x{t 2 )^^ . 


(15) 


Taking the averages by help of equation (lT4l) we arrive at 


-T-A 


AA(A) = 


(T- A)2 


dti 


x^(ti + A)^ 

L 

+ A)x{t2)‘j 


(16) 


With the new variable t' = t 2 — ti (assuming t 2 > ti) and by changing the order of 
integration we hnd the expression 


rT-A-r' 


V(A) = 


(T - A)= 


dr' dti (^x^{ti + A)^ 


x‘^{ti + r') 


(17) 


Now, the new variables x' = ti/A and y' = r'/A are introduced. Substituting equation 
([I]) into equation flTTD we obtain 


U{A) = 


167s:2A2«+2 


dy' 


rT/A-l-y' 


dx' 


(T-A)2 

X [{x’ + 1)2“ - 2{x' + lT{x' + vT + {x' + y'f^] . (18) 

Splitting the double integral over the variable x' into an integral over a square region 
and a triangular region yields 


dy' 


fT/A-2 


dx' + / dy' 


r^T/A-l-y' 


dx 


'T/A-2 


rTIA-2 


dx' / dy' + 


rT/A-l 


dx' 


fT/A-l-x' 


IT/A-2 


dy'. 


(19) 


From the double integrals from the power-law functions in equation flTSD . via equation 
(ITT)) we compute the nominator as 


16 A ' 2 A ^"+2 


(r — 1)2"+1 ^ (^ 3 q - -)- 1 )(^ 7 - _ X ) 2 o -|-2 


2r“+^(r - 1)“+^ 


2a+ 1 ' 2(a-t-l)2(2a + 1) 

r 2“+2 (^2a^ -f- a -f 1) 


(a+ 1)2 

/*r—1 


2(a + l)(2a + l) 2(a + l)2(2a + 1) 


dx' {x')'^~^^{x' + 1)^ 


0-1-1 


'0 


( 20 ) 
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in terms of the variable 
T 


T = 


A' 


The integral 


f*T—l 


h{T) = / dx' {x'+ ly 


( 21 ) 


( 22 ) 


remaining in the last term of this expression can, in principle, be represented in terms of 
the incomplete Beta-function. The denominator 'P(A) of the EB parameter ([H]) is just 
the squared time averaged MSD given by equation ([6]). We thus arrive at the expression 


P(A,r) = 




a+1 


(r - !)”+■) 


(23) 


L(a+l)(r-A) 

Note that the double analytical integration of equation (9) in [23] via Wolfram 
Mathematica yields a result, that is indistinguishable from equation 
demonstrated by the blue dots in £gure[3f3. 


as 


3.2. Expansions and Limiting Cases 

We here consider some limiting cases of the EB parameter based on expressions fl2UI) 
and fl20|) . In the limit a = 1 and for A/T <C 1 the leading order expansion in terms 
of A/T turns into equation (]9|). As it should the SBM process reduces to the ergodic 
behaviour of standard Brownian motion. 


3.2.1. The case 0 < a < 1/2. The general expression for the behaviour of the EB 
parameter in the range 0 < a < 1/2 follows from equation (12^ by help of the identity 
[equation (1.2.2.1) in Ref. 


x^{x -I- lydx = 


xP+\x + iy 


+ 


0. 


x^{x + iy ^dx, 


(24) 


p + q + 1 p + q + 1 
that can be checked by straight differentiation. Performing this sort of partial integration 
three times we reduce the power of the integrand so that in the limit r —?■ cxd the integral 
becomes a converging function. In the range 0<Q;<l/2we the End exact expression 


^ , . (r — l)“+^r“ Q!(r — l)“+^r“ ^ «(« — l)(r — l)"’''^r“ ^ 

“ 2(a + 1) ^ 2(a + l)(2a + l) ^ 4a(a + l)(2a + 1) 


a{a — 1)(q; — 2) 
4Q;(a -|- 1)(2 q; -|- 1) 


f*T — l 


X 


(x')“+n(r'+ 1)“-W. 


(25) 


The remaining converging integral can be represented in the limit A/T <C 1 via the 
Beta function: setting the upper integration limit (r — 1) —)■ cx) we obtain 


/ {xT^\x' + l)°-^dx' = B{a +2,1- 2a). 

Jo 

Then we arrive at the following scaling law for the EB parameter, 

^ 2a 

EB(«,A)~4C(«)(^ 


(26) 


(27) 
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Figure 3. Ergodicity breaking parameter EB of non-ageing SBM. (A) Results of 
numerical simulations are depicted by the data points. The analytical results based on 
equations (l20)) and ([23]) are given by the solid coloured lines. Data points for different 
lag times are shown in different colours. The values of EB for ultraslow SBM (IM]) at 
a = 0 and at a = Xj’l given by equation (|32|) are shown as the bigger black bullets, 
computed for A = 10°, 10^, and 10^. The larger orange bullets denote the same limits 
but without the additive constants to the leading functional dependencies with A/T. 
Parameters: the trace length is T = 10^, the number of traces used for averaging at 
each a value is N = 10°. (B) Exact and approximate analytical results for EB. The 
red, green, and blue curves are the exact evaluations of equation (|20)l . The dashed 
curve in the region a > 1/2 corresponds to equation (|3nil and the dashed curves for 
0 < Of < 1/2 are the results of [23]. The magenta curves in the region 0 < a < 1/2 
are according to the analytical expansion (EH) for given A values. The dark blue data 
points, coinciding with our exact result (1201) . follow from evaluating the double integral 
in equation (9) of [23] with Mathematica. 


where the coefficient is given by 

(1 — a)(2 — Q:)i?(a + 2,1 — 2 q:) — (2a^ + a + 1) 

“-2(a+l)y2a+l)-' 

The scaling form of EB versns (A/T) of eqnation (|27|) coincides with that proposed 
in reference [2l], and it is indeed valid for vanishing A/T and scaling exponents not 
too close to a = 0 and a = 1/2, see below. We hnd in addition that in the region 
0 < a < 1/2 the EB parameter of the SBM process becomes a sensitive fnnction of 
the lag time A, as shown in hgnre[3]A., both from onr theoretical resnlts and compnter 
simnlations. This means that no nniversal rescaled variable A/T exists, as is the case 
for standard Brownian motion. 

The asymptote fl27|) agrees with the resnlt (10) in [23] in the range 0 < a < 1/2 of 
the scaling exponent and for inhnitely large valnes r. Eqnation fl28ll above provides an 
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explicit form for the prefactor. In figure [Sfl the approximate expansion fl27p is shown 
as magenta curve. At realistic values A/T the asymptote fl27l) agrees neither with our 
exact expression (l20j) nor with the simulation data. As this demonstrates the exact 
expression (|20|) needs to be used a forteriori. The main reason is the finite r value used 
in the simulations: for very small A/T equation flTTll describes the exact result fl20ll 
signihcantly better (not shown). We note that away from the critical points at cc = 0 
and a = 1/2, equation fl27l) returns zero and inhnity, respectively (magenta curves in 
hgure[3B)- At these points special care is required when computing A in equation (|25l) . 
as discussed below. 


3.2.2. The case a > 1/2. For values a > 1/2 of the scaling exponent in the limit of 
small A/T the denominator fl2^ becomes 'T(r) ~ 4r^“. Note that here we need to 
include two more iterations of the integral in the last term of equation fl2^ by using 
equation (|2Tll . Then we arrive at a new integral term that is converging at r —?■ oo. 
Thus the nominator (|20D — after cancellation of the first three orders in the expansion 
in terms of large r—yields to leading order M{t) ~ 16a^r^"“^/[3(2Q; — 1)]. 

From the exact expression fl20l) by using the integration formula fl241) four times, 
we hnd the exact representation 


A(r) = 


(r — l)“'''^r“ Q:(r — l)“+^r“ ^ «(« — l)(r — l)"’''^r" ^ 


+ 


2(a + l) 2(a + l)(2a + l) 
a{a — l){a — 2)(r — i)"+ 27-"-3 
4a(a + l)(2a + 1)(2 q; — 1) 
a{a — 1)(q; — 2)(q; — 3) 


+ 


4a{a + l)(2a + 1) 


r*r—1 


, , y/ y/ "^dx. (29) 

4a(a +l)(2a +l)(2a - 1) Jo 

From this expression the leading term with the divergence at a = 1/2 is written 
explicitly and the remaining integral is converging only then. Plugging this expression 
into equations (Eup and fl23p and keeping terms of order ^ in the limit r S> 1 we 
recover the result of [23] given by equation fl30p . again valid in the range a > 1/2. 
Note that the divergence in the denominator of the last term in A in equation is 
compensated by the proper expansion of the remaining integral in Ji in the limit of large 
values of r for a > 1/2, see below. 

The EB parameter then scales as 

4 A 

lim EB(A) ~-. (30) 

This result coincides with expression (10) in [23] in the range a > 1/2. As mentioned 
already, special care is needed near the critical point 0 = 1/2. Equation (l30P implies that 
SBM is an ergodic process, with the EB parameter scaling strictly linearly with A/T 
as in relation (|9]) for Brownian motion, however, with an o—dependent prefactor of the 
form a^/{2a — l). In contrast to subdiffusive CTRW processes [3] SO] and heterogeneous 
diffusion processes izn the EB parameter for Brownian motion converges to zero and 
thus for sufficiently long measurement times the result of time averaged observables 
become reproducible. 
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3.2.3. The case a = 0 Now let us focus on the critical points a = 0 and a = 1/2 in 
detail. At a —)■ 0 the EB parameter of the ultraslow SBM process [36] can be obtained 
from equation (]20|) . To this end we hrst expand result (j20|) for small a using the identity 
_ goiog(a;) remaining integral Ji in equation (]22|) we hrst expand the integrand 


T-A 
0 


dt. 


in powers of small a and then integrate the expanded function in the limits 
The hrst two orders of the expansion in a in the nominator of EB disappear. Dividing 
the leading orders in in the nominator and denominator of EB and expanding for 
short lag times A/T <C 1 afterwards to the leading order we hnd 

4(7r^/6 — 1) 


lim EBusbm(A) 
A/T-i-O 


(log[T/A] + 1)^ 


(31) 


This result was obtained from independent considerations for ultraslow SBM as equation 
(20) in [36]. Note the logarithmic rather than the linear dependence of EB on A/T in 
this case, stemming from the ultraslow logarithmic scaling of the MSD and the time 
averaged MSD with (lag) time. 


3 . 2 . 4 . The case « = 1/2 Similarly, to explore the limit a —?■ 1/2 we hrst expand 
the exact result (1201) for A/'(A) in a around this point. In analogy to the case a = 0 
we expand the integrand in Ii in terms of powers of (a — 1/2) and then perform the 
integration over t from 0 to T — A. Dividing the expansion of the nominator fl20|) of 
EB, taken at a = i/i in the limit A/T —)■ 0, by the leading order of the denominator 
fl2^ in the same limit—scaling as 4r—we get 


lim EB„=i/2(A) 
A/T-s-O 


A 


log(T/A) + 2 log(2) - 5/6 


(32) 


The same expression can be obtained by expanding equation fl25l) valid in the region 
0 < a < 1/2. Alternatively result fj3^ can be obtained from the exact expression fl2^ 
valid for a > 1/2. In this case, however, due to a pole at a = 1/2 one more order in 
the power expansion near 0 = 1/2 needs to be properly evaluated when expanding A. 
Then, the divergence in the denominator of the prefactor of the last term in equation 
fl2^ becomes eliminated and the EB parameter stays continuous as o —)■ 1/2. 

Compared to the case o = 1 of Brownian motion the result 0321) for EB features 
a weak logarithmic dependence on A/T. As expected the values of EB according to 
equation (j32l) are very close to the exact solution (j20|) . as shown by the larger black 
bullets for a = 1/2 in £gure[3]A. Note that for finite T/A values the additional constants 
following the leading functional dependencies in equation fl3Tl) and equation fl3^ play a 
signihcant role, as seen in hgure [3]A. The agreement of these EB values with the exact 
predictions of equation fl20l) and computer simulations is particularly good for smaller 
A/T values, as expected based on the large r expansions used in the derivation of 
equations (l3TD and (1321) . 


II With regard to the higher order expansion taken below, this corresponds formally to an expansion 
of order {a — 1/2)°. 
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3.3. Computer Simulations 

We implement the same algorithms for the iterative computation of the particle 
displacement x(t) as developed for the heterogeneous diffusion process [2T] and the 
combined heterogeneous diffusion-scaled Brownian motion process [26]. We simulate 
the one dimensional overdamped Langevin equation 

'SSL = VW X «*) (33) 

driven by the Gaussian white noise f{t) of unit intensity and zero mean. At step i + 1 
the particle displacement is 

Xi+i -Xi = ^/2[D{ti) + C]{yi+i - yf), (34) 

where the increments (|/i+i — |/i) of the Wiener process represent a 5 correlated Gaussian 
noise with unit variance and zero mean. Unit time intervals separate consecutive 
iteration steps. To avoid a possible particle trapping at the pole of D{t) we introduced 
the small constant C = 10“^ in analogy to the procedure for heterogeneous diffusion 
processes [2I|. The initial position of the particle is Xq = x{t = 0) = 0.1. 

Our simulations results shown in hgure [3]A conhrm the validity of the general 
analytical expressions (1201) and fl2^ making up the EB parameter in the whole range 
of the scaling exponent a. We also find that the short lag time expansion fl30|) agrees 
well with the exact solution and simulations at a > 1/2 (hgure EB)- In the range 
a > 1/2 the EB parameter for A/T -C 1 is nearly insensitive to the lag time and grows 
with a in accord with equation fl30|) . In particular, the full analytical expression for 
EB (equations (12(11) and fl2^ ) and the results of the simulations show no divergence at 
a = 1/2, in contrast to the approximate results of reference [24] . 

Figure EB- also shows the approximate EB values (j3Tl) for ultraslow SBM as well 
as EB at a = 1/2 from equation (l32|) indicated as larger points. These points are close 
to our predictions for SBM at a —)■ 0, in particular, for small A/T values when the 
approximations used in deriving the corresponding equations are better satished. As 
the ratio A/T grows and the scaling exponent converges to zero, a 0—indicating 

progressively slower diffusion—the results of our simulations start to deviate from the 
exact analytical results (l20|) and fl23l) . as shown in hgure El In this limit apparently 
better statistics are needed in the simulations. 

In hgure m we show that EB scales with the trace length T approximately as 1/T^“ 
for 0 < a < 1/2 and as 1/T for a > 1/2; compare to the results in hgure 1 of reference 

4. Ergodicity breaking of ageing scaled Brownian motion 

We consider the ergodic properties of ageing SBM, where ta denotes the time span in 
between the initiation of the system and start of the measurement. The ergodicity 
breaking parameter is dehned through the ageing time averaged MSD (compare 
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Figure 4. EB parameter for non-ageing SBM versus trace length T. The solid lines 
represent the exact results according to equation (1201) . Parameters: A = 10 and 
N = lO^. 


equations ffTOj) and flTT]) ) as 


EB.(A) 


(sK^f) - (^;(A))' 
(f(A))' 


A/'a(A, r) 
I^a(A,r) 


For the numerator we find in full analogy to the non-ageing situation 

X + A)) - (x(ti + A)i((2))]^ . 


(35) 


(36) 


Changing the variables as above for the non-ageing scenario, r' = t 2 —ti, we switch the 
limits of integration using ti(r') = T -|- — A — r' and then split the integrals over r' 

to compute the pair correlators using the property flT^ . This yields the representation 
of the nominator of EB in terms of one-point averages only, 


A4(A) = 


(T - A)2 

rA pT+ta-A-r' 


X 


dr' 


dti 


'ta 


^ 2 


{x^iti + A)) - {x^iti + t')) . (37) 


We proceed by inserting the MSDs of equation ([T]) and arrive at 




^T/A+ta/A-l-y' 


dx' 


{x' + 1)2“ - 2(x' + lT{x' + y'T + {x' + y') 


/\2a 


X 


(38) 
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Changing the order of integration and splitting the integral over x' we get in terms of 
the variables t = T/A and 
ta 


Ta = 


A 


(39) 


that 


A/'a(A, r) = 


16iC2A2“+2 
(T - A)2 j 

+ (x' +!/')'“ 


-\-Ta — 2 


dx' / dy' {x' + ly°'— 2{x' + lYi^x'+ y'Y 


rr+Ta — l 


+ 


dx' 


I'T+Ta—l — x' 


dy' 


' T+Ta -2 


X 


(x' + 1)’“ - 2(x’ + l)“(a:' + j/')“ + (i' + y'f^ 


(40) 


Finally, taking the integrals in the nominator of EB for ageing SBM yields 


A/'a(A,r) = 


16iC2A2«+2 

(T - A)2 


{r + Tg- 1) 

2Q! + 1 


20+1 


('^a + 1) 


2q+1 


2q; + 1 


(3a + l)(r + Ta - l)2"+2 ^ (3 q; + + 1 ^ 01+2 


2(2a + l)(a + l)2 


2(2a + l)(a + l)^ 


2(2a + l)(a + 1) 

('r + ra)^“+^ 2(r + ra)“+^(r + Ta - 1)“+^ 


2(2a + l)(a + 1) 


(a+ 1)2 


OT+Ta-l 


dx' {x'T+\x' + iy 


-f-1 


(41) 


Here we again denote 

h{T,Tg) = 


r-T+Ta-l 


dx'{x')‘^^\x' + 1Y 


(42) 


The denominator of EB follows from the time averaged MSD flTT]) . namely [26[ US] 


Vg{A,T)= (SKA)) = 


2iF„A"+^ 


(r + - {Ta + 1) 


a+1 


(a + l)(T-A) L' 

-(r + T.-l)“+‘+T“+‘])'. (43) 

The final EB breaking parameter fl5^ for ageing SBM turns into expression fl2U]) for the 
non-ageing case, Tg = 0. 

In the limit of strong ageing, Tg ^ T ^ A, the time averaged MSD scales as 

(44) 


2aKXg~ A 


and the nominator of EB grows as 


A4(A,r) ~ 16A'>=“T-"(+Tf-+/3) 


(45) 


to leading order in large Tg values and long trajectories. Then, the ergodicity breaking 
parameter follows the Brownian law ([9]). This limiting behaviour is supported by the 
simulations of strongly ageing SBM shown in hgure [5l Moreover, it is similar to that of 
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a-1 


Figure 5. EB parameter for ageing SBM. Results of simulations are shown by 
the points and the analytical results dm are represented by the solid lines of the 
corresponding colour. Parameters: A = 10, T = 10"^, and N = 10^. 


ageing ultraslow SBM [36] . Physically, in the limit of long ageing times Tq the diffusivity 
D{t) changes only marginally on the time scale T t of the particle diffusion, so that 
the entire process stays approximately ergodic. 

In the opposite limit of weak ageing, Tq T, we observe that {^5^ (A, r)) ~ 
+ aTaT°'~‘^), and the nominator of EB to leading order of short Ta and 
long T values produces A/'(A,r) ~ 16iP^A^“r“^(Q;^r^"“^/[3(2a — 1)]). Consequently 
the EB parameter to leading order is independent of the ageing time and follows 
equation fl30l) as long as a > 1/2. 

Figure [5] shows the simulations results based on the stochastic Langevin process of 
ageing SBM. We hnd that in the limit of strong ageing, consistent with our theoretical 
results the EB of ageing SBM indeed approaches the Brownian limit ([9|). For weak 
and intermediate ageing the general EB expression (IM]) is in good agreement with the 
simulations results, compare the data sets in hgure [5l Finally hgure |6] depicts the 
graph of EB versus ageing time explicitly, together with the theoretical results fimj) 
and (|43|) . We observe that EB decreases with the ageing time and this reduction is 
particularly pronounced for strongly subdiffusive SBM processes. The latter also feature 
some instabilities upon the numerical solution of the stochastic equation for long ageing 
times. 
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CD 

LU 



Figure 6. EB parameter for ageing SBM versus ageing time ta- Analytical results (1411) 
and (l43l) for different a values are represented by the solid lines. Some instabilities in 
the simulations are visible at long ageing times, in particular for small a. Parameters: 
A = 10, r = 10^, and N = lO^. 


5. Conclusions 

We here studied in detail the ergodic properties of SBM with its power-law time 
dependent diffusivity D{t) ~ In particular, we derived the higher order time 

averaged moments and obtained the ergodicity breaking parameter of SBM, which 
quantihes the degree of irreproducibility of time averaged observables of a stochastic 
process. For the highly non-stationary, out-of-equilibrium SBM process we analysed the 
EB parameter with respect to the scaling exponent a, the lag time A, and the trace 
length T. We revealed a non-monotonic dependence EB(a). In particular, we showed 
that there is no divergence at a = 1/2, in contrast to the approximate results of |24j . 
We also obtained a peculiar dependence for the EB dependence on the trace length T, 
EB(T) ~ 1/T^“ for 0 < a < 1/2 and EB(T) ~ 1/T for a > 1/2, in agreement with jMj. 
We also obtained analytical and numerical results for EB for ageing SBM as function of 
the model parameters and the ageing time ta- 

Our exact analytical results are fully supported by stochastic simulations. We hnd 
that over the range a > 1/2 and for A/T 1 the EB dependence on the lag time and 
trace length involves the universal variable 1/r = A/T, as witnessed by equation fl30|) . 
For arbitrary lag times and trace lengths the general result for ageing and non-ageing 
SBM are, however, more complex, see equations (l20l) and 6a. These are the main 
results of the current work. For strongly subdiffusive SBM in the range of exponents 
0 < a < 1/2 the ergodic properties are, in contrast, strongly dependent on the lag time 
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A. The correct limit of our exact result (I2U]) was obtained for the EB parameter of 
ultraslow SBM with a —)■ 0 and for SBM with exponent a = 1/2. Although EB has 
some additional logarithmic scaling at this point, it reveals no divergence as a = 1/2 is 
approached. 

We are conhdent that the strategies for obtaining higher order time averaged 
moments developed herein will be useful for the analysis of other anomalous diffusion 
processes, in particular for the analysis of hnite time corrections of EB for fractional 
Brownian motion 051 or for processes with spatially and temporally random diffusivities 

[501 ED. 
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